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Abstract 

We consider Af = 3 supersymmetric Chern-Simons (CS) theories that contain prod- 
uct U(N) gauge groups and bifundamental matter fields. Using the matrix model of 
Kapustin, Willett and Yaakov, we examine the Euclidean partition function of these 
theories on an S 3 in the large N limit. We show that the only such CS theories for 
which the long range forces between the eigenvalues cancel have quivers which are in 
one-to-one correspondence with the simply laced affine Dynkin diagrams. As the A n 
series was studied in detail before, in this paper we compute the partition function for 
the D4 quiver. The D4 example gives further evidence for a conjecture that the saddle 
point eigenvalue distribution is determined by the distribution of gauge invariant chi- 
ral operators. We also see that the partition function is invariant under a generalized 
Seiberg duality for CS theories. 



Contents 



1 Introduction 



1 



2 TV = 3 SUSY CS Theories and ADE Dynkin diagrams 



4 



3 Generalized Seiberg Duality 



8 



4 Z?4 type quiver 

4.1 Comparison with the Numerical Saddle Point 

4.2 Operator counting 



m 

M 



5 Discussion 



HI 



A A Simple Class of Solutions 



El 



1 Introduction 

We hope that supersymmetric (SUSY) gauge theories in 2+1 dimensions will help us learn 
about general features of 2+1 dimensional gauge theories which in turn might shed light on 
certain condensed matter systems with emergent gauge symmetry at low temperatures. For 
many years, there has been strong evidence that, similar to the four dimensional case, three 
dimensional gauge theories are related to each other by a large number of non-perturbative 
dualities. In four dimensions, anomalies provided quantitative tests of these conjectured 
dualities, while in three dimensions, because of the absence of anomalies, the situation was 
more difficult. Recently, with the introduction of methods for calculating a superconformal 
index on S 2 xS 1 [1] and the Euclidean partition function on5 3 [2] , the situation has changed 
dramatically. 

In this paper, we continue an investigation, started in [31 13 [5], of the the large N limit 
of the S" 3 partition function of supersymmetric (SUSY) Chern-Simons (CS) theories. The 
partition function Z S 3 is calculated using the matrix model derived in ref. [2] by localization 
(later improved by [H [7] to allow matter fields to acquire anomalous dimensions) . For the 
CS theory at its superconformal fixed point, the matrix model of [2 J computes exactly 
the partition function and certain supersymmetric Wilson loop expectation values. The 
theories we examine here have M = 3 SUSY, a product U(N) gauge group structure, and 
field content summarized by a quiver diagram. While in our previous work [31 HI [5] we 
examined CS theories that had a U(N) d gauge group, in this work we relax the constraint 
that the ranks of the U(N) groups all be equal to each other. 

We examine quantitatively some of these three dimensional dualities. We will study the 
generalized Seiberg duality of [HIE]- Modulo convergence issues, refs. PUHH] demonstrated 
that \Zgs\ is invariant under this duality. In this paper, we investigate what this invariance 
looks like in the large N limit. We will also examine the AdS/CFT duality between these 
gauge theories and Freund- Rubin compactifications of eleven dimensional supergravity. 

A big motivation for this paper and our previous work [31 HI [5] is the hope that the 
matrix model can shed light on the microscopic origin of the mysterious iV 3 / 2 scaling of the 
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free energy |12] indicated by AdS/CFT. An early investigation of this matrix model j!3| 
demonstrated that the free energy of maximally supersymmetric SU(N) Yang-Mills theory 



While one could declare victory, the localization procedure of [2] leaves the microscopic 
physics of this scaling obscure. Here we study a larger class of gauge theories using different 
methods in the hopes of shedding further light on this large N scaling^] 

In the context of the AdS/CFT correspondence, the free energy F can be computed 
from a classical gravity model dual to the strongly interacting, large N field theory. For a 
CFT dual to AdS^ of radius L and effective four-dimensional Newton constant Gn, F is 
given by [15] 



These types of AdS^ backgrounds arise as Freund-Rubin compactifications of AdS± x Y 
of M-theory, where Y is a seven-dimensional Sasaki-Einstein space supported by N units 
of four-form flux. The quantization of L in Planck units implies that at large iV eq. ^ 
becomes [3] 



(Here, the volume of Y is computed with an Einstein metric that satisfies the normalization 
condition R mn = 6g mn -) These types of Freund-Rubin solutions arise as the near horizon 
limit of a stack of iV M2-branes placed at the tip of the Calabi-Yau cone X over Y. The 
CFT in question is then the low energy field theory description of the M2-branes. The 
matrix model can be used to compute not only the three-halves scaling but also the volume 
factor and thus allows for quantitative tests of proposed AdS/CFT dual pairs. 

Note that the conjecture of ref. [12] preceded the discovery of AdS/CFT and also con- 
cerned the scaling not of F but of the thermal free energy. AdS / CFT makes clear that the 
iV 3 / 2 scaling comes from the over-all normalization of the gravitational action and thus is 
a universal property of many quantities in the dual field theory. 

In trying to elucidate the connection between matrix model quantities and field theory 
quantities in [H [3] , we found an interesting relationship between the saddle point eigenvalue 
distribution and the distribution of gauge invariant operators in the chiral ring. We now 
review this relationship. The localization procedure of [2] reduces the partition function 
to an integral over d constant N x N matrices a a where a a is the real scalar that belongs 
to the same M = 2 multiplet as the gauge connection. Ref. [3] presented a procedure 
(later improved in [16]) to evaluate the matrix model by saddle point integration in the 
large N limit. At the saddle point, the real parts of the eigenvalues A^- of a a grow as a 
positive power of N while their imaginary parts stay of order one as iV is taken to infinity. 
Additionally, the real parts of the eigenvalues are the same for each gauge group. Therefore, 

1 See [14] for a recent Fermi gas interpretation of the microscopic physics. 





(3) 
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in order to find the saddle point, one can consider the large iV ansatz 

\f = N a Xj + iy ad + ... (4) 
As one takes N — > oo, the xj and y a j become dense; one can introduce an eigenvalue density 



1 N 

p W = jS&on'E 5 ^-^ ■ (5) 



.r- 



The chiral ring of these quiver gauge theories consists of gauge invariant products of the 
bifundamental and fundamental matter fields and monopole operators modulo superpoten- 
tial and monopole relations. At large N, only the so-called "diagonal monopole operators" 
are important. These operators turn on the same number of units of flux through the di- 
agonal U(l) subgroup of each U{N) gauge group. Operators in the chiral ring therefore 
have an associated R-charge r and a (diagonal) monopole charge m. To characterize these 
operators, we introduce the function tjj(r, m) which counts the number of operators in the 
chiral ring with R-charge less than r and monopole charge less than to. In [H [5], we found 
evidence for the following relation between the saddle point eigenvalue distribution and 
il){r, to): 

d 3 ip 



dr 2 d' 



m 



= r ~p(x) , (6) 

m=rx/n r 1 

where \x is a Lagrange multiplier that enforces the constraint J p(x)dx = 1. Defining a 
function tpx( r , r n) that counts the number of operators without the field X, there is a 
similar conjectured relation [HE] between ifi^ i r i to) and the y a j which we shall review in 
section 14.21 

In the current paper, we begin in section [2] by examining CS theories with N = 3 
SUSY with matter fields only in bifundamental representations of the product gauge group. 
However, we relax the constraint that the ranks of the U(N a ) gauge groups all be equal. 
We find that the saddle point procedure described in [3] for evaluating the matrix model 
works if the ranks of the gauge groups satisfy the following condition: 

2N a = ]T N b , (7) 
b|(o,6)eJ5 

where E is the set of bifundamental fields. (In other words, the number of flavors for a 
given gauge group should be twice the number of colors.) If this condition is not satisfied, 
there will in general be long range forces between the eigenvalues that render the method 
in [3] inapplicable. Interestingly, this condition implies that the gauge groups and matter 
content are described by a simply laced, affine Dynkin diagram (see figure [I]). 

For a 3 + 1 dimensional theory, the condition (J7| would seem very natural because it 
implies the vanishing of the NSVZ beta function and the existence of an IR superconformal 
fixed point. In this 2 + 1 dimensional context, it is less clear where such a condition comes 
from. Gaiotto and Witten [T7] call such quivers balanced. In their study of M = 4 Yang- 
Mills theories in 2+1 dimensions (without CS terms), they find that the Coulomb branch 
of the moduli space for balanced quivers has an enhanced global symmetry. 
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Figure 1: The simply laced affine Dynkin diagrams of ADE type. 

When the ranks N a of the gauge groups are all equal, the Dynkin diagram is of A n 
type and was studied in detail in our earlier papers IU [5] . However, one may also try to 
solve the matrix model for the D n and E n quivers. In section [3j we discuss how generalized 
Seiberg duality acts on the A n , D n , and E n quivers in the large iV limit. In particular we 
investigate how the saddle point eigenvalue distribution transforms under Seiberg duality. 
In section |4j we investigate the D4 case in detail. We find that the relation ^ between p(x) 
and the chiral ring continues to hold. A discussion section contains some partial results and 
conjectures. In appendix [Aj we find a simple class of solutions for any ADE quiver. 

Note Added: After this paper appeared, we became aware of [18] . Section 4 of [T8] has 
some overlap with section [3] of this paper. 

2 M = 3 SUSY CS Theories and ADE Dynkin diagrams 

We are interested in Chern-Simons theories with matter and at least J\f = 3 supersymmetry 
in 2+1 dimensions. We assume the gauge group to have the product structure Ylt=i U(N a ). 
In other words, for each U(N a ) factor, there exists a vector multiplet transforming in the 
corresponding adjoint representation. We denote the Chern-Simons level for each U(N a ) 
factor k a . In addition to the vector multiplets, we allow the theory to contain hypermul- 
tiplets in bifundamental representations of the gauge groups. More precisely, in N = 2 
language, the hypermultiplet consists of a pair of chiral multiplets (X a b, Xb a ) transforming 
in conjugate bifundamental representations, (N a ,N{,) and (Nb,N a ). The field content of 
such a theory is conveniently described with a quiver. For each U(N a ), we draw a vertex. 
For each pair (X a b, Afe a ), we join nodes a and b with an edge. Let E be the set of edges. 
Recall that an M = 3 vector multiplet contains an M = 2 vector multiplet V and a 
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j\f = 2 chiral multiplet <I> transforming in the adjoint representation of the gauge group. 
The scalar component (j) of $ combines with the auxiliary scalar a in V to form a triplet 
under the SU(2)r symmetry. Similarly, the pair (X a 6, X& a ) form a doublet. The j\f = 2 
formalism leaves only a J7(1).r C SU(2)r of the R-symmetry manifest. Under this U(1)r, 
$ has charge one and (A^A^) has charge 1/2. 

We would like to compute the Euclidean partition function of such an j\f = 3 Chern- 
Simons theory on an S 3 . As explained by Kapustin, Willett, and Yaakov [2j, this path 
integral localizes to configurations where the scalars a a in the vector multiplets are constant 
Hermitian matrices. Denoting the eigenvalues of a a by A aj j, 1 < i < N a , the partition 
function takes the form of the eigenvalue integral 



/ l]JdX a A L v ({X a>i })L m ({X a>i }) = j l]JdX a , 

\ a,i J \ a,i 



Z = J \\l dX ^ I L v{{KA) L m{{K,\) -- I ( | | '/A,,., I cx|) [-./-'({A,,.,})] . !8) 
where the vector multiplets contribute 

L - = II J^l ( II 2 8inh K A <M - Aa,,)] J exp ( iTT ^ fc a A^- J (9) 

a=l a ' \i>i J \ a,j J 

and the bifundamental matter fields contribute 



Lm ~ n n 2cosh[7r(A ^_ A )i • ( io ) 

(a,b)eE i,J 



We wish to approximate this integral using the saddle point method in the limit where 
all of the N a are large. In particular, we will assume that N a = n a N where the n a are 
relatively prime integers and iV S> 1. The saddle point equation for the eigenvalue A 0) j is 

^ — = 2ik a Xg t i + y~] 2 coth[7r(A a) j - X a j)] - V] V] tanh[7r(A aii - A 6j -)] = . (11) 

i/i 6|(a,6)eB J 

Guided by numerical experiments, we will assume that we can pass to a continuum 
description where each set of N a eigenvalues lie along n a curves, iV eigenvalues per curve: 

X a ,i(x) = N a x + iy aJ (x) , (12) 

where a > and I = 1, . . . , n a . The density of eigenvalues p(x) is assumed to be independent 
of a and /, and we normalize the density such that J p(x) dx = 1. At large N, we may use 
the approximations 

coth[7r(A a!i - Xb,j)] ~ tanh[7r(A aj i - X b j)] w sgn(xj - sc 3 -) . (13) 

Assembling these large N approximations, the leading order term in the saddle point 



equation (11) is (assuming a < 1) 

2n a 



rib \ N dx' sgn(rc — x')p(x') . (14) 
b\(a,b)eE j 
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For this term to vanish, a sufficient condition is that 



2n a = ^ n b . 

b\(a,b)EE 



(15) 



To find saddle points for more general theories, one obvious thing to try is a more general 
ansatz for the eigenvalue distribution, in particular one that does not assume p(x) is the 
same for each gauge group. The condition (15) is surprisingly stringent. The only allowed 
quivers that satisfy this condition are the extended ADE Dynkin diagrams! In this context, 
the n a are called the marks. To each node of the Dynkin diagram, we associate a simple 
root a a . For an affine Dynkin diagram, we have one additional root 9 associated to the 
extra node. This extra root, which is the highest root in the adjoint representation, satisfies 
the condition 9 = Y^ a n a a a . The mark associated with 9 is defined to be equal to one. In 
the rest of the paper we will assume that (15) holds. 



Next we examine the large iV limit of F itself, assuming ( 15 ). It is convenient to separate 

(16) 



the free energy into pieces: F = F tree + -^ioop + ln(iV a !) where 

F tre e = -»7T ^ k a \ 2 a i . 



a. i 



The X^a m (-^') contribution will turn out to be subleading in N, and we hereby drop it. 
In the large N limit, we find that 



tree 



o=l 



-iN 2a n, 



/W-a /» 
x 2 p{x)dx + 2N a K j xy a j{x)p(x)dx + 0{l) 
i=i ^ 



■ (17) 



The loop contribution to F is 



F\ 



loop 



N 2 j dxdx' p(x)p(x') x 

/ d n a n a 



^2 ^2 hi2 sinh[7riV a |x - x'\ + m{y a>I {x) - y a> j(x')) sgn(x - x')] 



a=l 1=1 J=l 

n a n b 



+ Yl ^^ln2cosh[7rAr a (x- 2 ; , ) + ^(ya,/(3 ; ) -yb,j{x'))\ • 

(a,b)£E I=1J=1 J 



(18) 



This integral is dominated by the region x ~ x' which suggests we change variables from 
x' to £ = 7rN a (x — x') so that to a first approximation at large N we may write y a (x') = 
y a (x — N^ a ^/ir) fa y a (x). The integral reduces to 



Fi 



N 



2-a 



oop 



7T 



dx d£ p(x) 2 x 



d n a n a 

~^2^2^2 ln2 Sillh [1^1 + * ^(^(Vajix) - Va,j(x))] 
a=l 1=1 J=l 

n a n b 

(a,b)eE 1=1 J=l 



(19) 



6 



Given (15), the integral over £ can be performed and is convergent. The following two 
In 2 sinh [\x\ + sgn(x)ia] dx — M 2 



intermediate results are useful 



lim 

M-5-oo 



M 

lim 

M-^oo 



I In 2 cosh(x + id] dx - M 2 

J-M 



1 

12 ~ 4 

Vr2 _ 1 

12 4 



arg(e 2 -) 2 



The end result is that 



N 



2-a 



oop 



47T 



dx p(x) 



d n a n a 

E £ £ arg i ' 

,a=l 7=1 J=l 

- E EE a "i' 

(a,6)e-B Z=l J=l 



^i(y a t-Vb ,j) 



(20) 



To find a nontrivial saddle point, we have two choices. We can balance N 01 against 



the first term in (17), in which case a = 1/3. Alternately, we can assume that Yl a na ^ a = ^ 

171) in which case a = 1/2. The first case 
We shall make the second choice, leading 



and balance N 2 ~ a against the second term of 
has a massive type IIA supergravity dual [16 
to an eleven dimensional supergravity dual. 

The free energy that we must maximize is then 



F 



N 3 ' 2 / p(x) 



a 1=1 



2vrx e E + ^ ( E E E ar § 



are e 



va=l 7=1 J=l 



,2ni(Va,I-Va,. 7-1/2) 



E EE arff!r 

(o,6)eB/=l J=l 



^-^(J^-l) (21) 



where we have added a Lagrange multipler /x to enforce the constraint that the eigenvalue 
density integrates to one. 

We believe that the free energy can always be written in a way such that it depends 
only on the average values 

5, = -Ew (22) 
n a — 

rather than on the y a j individually. In particular, we conjecture that F will take the form 



F = ttN 3/2 I p(x) 2x E k a n a y a + p(x) E n ^ \\ ~ (Sa ~ Vbf J 

J n f~ U\r- TP \ / 



-2^ (/*,)*-! 



(a,b)£E 



dx 



(23) 



We have verified this conjecture in the case of the A n and D4 quivers. (Indeed, in the A n 
case, this result was given in [HE].) Note we are not claiming that the differences y a j — yb,j 



2 The ansatz ( |l2| needs to be slightly modified in this case to allow both the real and imaginary parts of 
A to scale with JV 5 ^ 3 . See \M- 
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remain bounded between ±1/2. Indeed, for the D4 case, the differences are not bounded. 
We are arguing that a more subtle cancellation happens among the integer shifts introduced 
in evaluating the arg functions such that (21) reduces to (|23[)p| 



3 Generalized Seiberg Duality 

Modulo some subtle convergence issues, refs. [TOl [11] demonstrated that the absolute value 
I Z s s I of the partition function Q is invariant under a generalized Seiberg duality [U [9] . 
For a U(N C ) theory with Nf flavors and CS level k, this generalized Seiberg duality acts by 
sending N c —> Nf — N c + \k\ and k — > —k. In the context of our quiver theories, Seiberg 
duality can be performed on any individual node of the quiver. The bifundamental fields 
look like flavors from the point of view of the node to be dualized. We have the long range 



force cancellation condition (15) which implies that Nf = 2N C . We also work in the large 
N limit which implies that k <C Nf,N c . Thus, under Seiberg duality N c is invariant to 
leading order in N. 

In the quiver context, in addition to changing the CS level of the dualized node, the 
duality also changes the CS levels of the neighboring nodes. If we dualize node a, then 
for b\(a,b) £ E, kt, — >■ + k a . Note that the sum ^2 a n a k a is invariant under the duality 



given (15). To motivate why the CS levels change in the way they do, we recall the brane 
construction of the A n quivers |20t |2"T] . We start with type IIB string theory with the 6 
direction periodically identified. We place N D3-branes in the 0126 directions. Intersecting 
these D3-branes at intervals in the 6 direction, we insert (l,p a ) 5-branes that span the 012 
directions as well as lines in the 37, 48, and 59 planes that make angles 9 a = arg(l + ip a ) 
with the 3, 4, and 5, axes respectively. The choice of 9 a guarantees that the construction has 
the right number of supercharges for the 2+1 dimensional gauge theory living on the D3- 
branes to have N = 3 supersymmetry. Each unbroken D3-brane interval in the 6 direction 
corresponds to a U(N) gauge group in the A n quiver. The CS level of the gauge group is 
given by the difference p a — p a +\. Seiberg duality corresponds to interchanging neighboring 
(l,p) 5-branes j^] 

For the simply laced affine Dynkin diagrams A n , D n , and E n , Seiberg duality can be 
reinterpreted as the action of the Weyl group. (There is a similar story for 3+1 dimensional 
gauge theories with M = 1 SUSY [22].) To n of the gauge groups we associate the corre- 
sponding simple root a a £ W 1 of the Dynkin diagram. The n + 1st gauge group is assigned 
the root a n +i = —9 where 9 is the highest root. The CS level is k a = a a ■ p where p is an 
arbitrary weight vector. Note that 

n 

} J n a k a = -9 ■ p + S ^n a a a ■ p = , (24) 

a a=l 

since Yla=i n a a a = 9. (The mark associated with 9 is defined to be one.) 



3 The second term of «2lb and the second term of (23 1 do not have the same periodicity properties as 



functions of the y a ,i- In the limit y a — y b — > 00, we do not expect (21l and (231 to be equal. What we 



are claiming is a region in y a j space where (21 1 and (231 are equal and which contains the saddle point 
eigenvalue distribution. 

4 A similar brane construction involving an orbifold should exist for the D n quivers [T71 . 
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Returning to the brane construction for the A n quivers, the Weyl group is the permu- 
tation group acting on the (l,p a ) branes. We start by assigning to the n + 1 gauge groups 
the n simple roots a a = e a — e a +i with e a a unit coordinate vector in IR n+1 and the highest 
root —9 = e n +i — e%. Note that these n + 1 roots live in a M. n subspace of M n+1 . Let 
V = (pi)P2> • • • ;Pn+i) be a weight vector corresponding to the charges of the 5-branes. As 
claimed, the CS levels are then k a = a a ■ p. 

The results of refs. [10} [TT] imply that \Zg3\ is invariant under this generalization of 



Seiberg duality for any N. The large N limit (21) must also be invariant, but it is still 



interesting to see how the invariance arises in our case. If we shift k a — )• —k a , kb —> kb + k a 



for b\(a, b) G E, then the first term of (21) can be made invariant if we take 

n b 

v*,i - y^ 1 + Yl Yl yb > j ■ ( 25 ) 

/ / b\(a,b)EE J=l 



In the case when n a = 1, it is clear that not just the first term of (21) but all of (21) is 
invariant under such a transformation. For example, for the A n quivers, the transformation 
on y a j reduces to 

Va -> -y a + y a -i+ y a +i ■ (26) 

To understand Seiberg duality on the internal nodes of the D n and E n quivers, we should 
describe how each y a j transforms rather than how the average transforms. However, it is 



straightforward to demonstrate that the conjectured form of the free energy ( 23 ) is invariant 



under (25). Thus, if we believe (23), i.e. that the free energy depends only on the average 
values y a , then F is invariant under our generalized Seiberg duality. 

We should note that our large- iV approximation is insensitive to shifting all of the y's 



by the same amount. The true Seiberg duality is actually a combination of (25) with such 



a shift. (See section 4.1 for an example of this phenomenon. 



4 D4 type quiver 

Let us examine the D4 type quiver. Let the external vertices have CS levels k a , a = 1, . . . , 4 
and the internal vertex have CS level k where Ik = — Y2 a k a - We label the imaginary part of 
the eigenvalues of the internal vertex yo^x) and 2/0,2(2;) and define 5y a j(x) = y a (x)—yoj(x). 
Based on numerical simulation and intuition from the A- type quivers [5J, we expect the 
saddle point eigenvalue distributions are piecewise smooth functions of x. 

The action of generalized Seiberg duality separates the space of CS levels into different 
Weyl chambers. It is enough to analyze the theory in a given Weyl chamber. Instead we 
will make the choice k\ > k2 > k^ > k^> Q which includes three different Weyl chambers. 
At the end, we will be able to see explicitly how Seiberg duality acts on the solution. 

In the central region \x\ < n/2k\, the equations of motion yield 

,W = f; W.) = |. (27) 

Once \x\ > fi/2\k\\, the solution gets more complicated. In general one of the #2/1, /> which 
for convenience we choose to be Syi t %, will saturate at 1/2 and the other will continue to 
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grow such that Syi i2 > 1/2. There are three sub-cases to consider: 1) k\ — k 4 > k 2 + k 3 , 2) 
h + k 4 > k 2 + k 3 > ki - k 4 and 3) k 2 + k 3 > k\ + k 4 . Let 5y a = y a - (y ,i + 2/0,2)/2. The 
next three piecewise linear regions of the solution are 

2fci _ _ ifei + & 2 2 ' v ; 

1 r fc 2 x e k 3 x k 4 x 

2/o,i - 2/0,2 = — 1 , 



k 4 + k 2 ki + k 3 



1 1 fc 3 X fc 4 X 

<W = 2, «W,2=2. to = ^, ^4 = ^ 
(fci - k 2 )x 



2/0,1 - 2/o,2 



2p 



" <x<x 4 : p= ^_ (2fci + fa + fcs)x > (30) 



^i + fe 3 ~~ r 2 2 

02/i,i = 2 ' ^^2.2 = 2 ' y3 ' 2 = 2 ' V4 = ~2p ' 

1 , (h- k 2 ~ h)x 
2/0,1-2/0,2=3 + - , 

where in cases 1 and 2, x 4 = n/{k\ + k 4 ) and in case 3, x 4 = n/(k 2 + k 3 ). The fifth and 
final piecewise linear region of the solution is different in each of the three cases. In cases 
1 and 2, we find 

M ^ ^ o (3&i + /c 2 + /c 3 + fc 4 )x 

< x < x 5 : p = 2/z , (31) 



ki + k 4 ~ ~ 2 

<fyl,l = 2 ' ^2,2 = ^ > ^2/3,2 = 2 ' <%4,2 = " , 

1 (fe + fc 2 + &3 + A; 4 )x 
2/o,i - 2/0,2 = 2 ^ , 

where X5 = /i//ci in case 1 and yo,i — 2/0,2 — > 1 at the endpoint. In constrast, in case 2, 
X5 = 2fi/ (k\ + k 2 + k 3 + k 4 ) and yo,i — 2/0,2 — ^ at the endpoint. In case 3, we have instead 

<x<- — ^- — : p = 2/x - (fei + fc 2 + h)x , (32) 



k 2 + k 3 k 4 + k 2 + k 3 + k 4 

~ 1 c 1 _ 1 r k^-X 

oyi,i = - , 02/2,2 = - , dy 3:2 = - , dy 4il = — , 

2/o,i - 2/0,2 = . 

In the three cases, the value of the Lagrange multiplier is 

2 6 2k 1 1 1 

case 1 : —5 = -. — V To ~ 1 ; 1 ; 1 r~ > ( 33 ) 

/i^ Ki fcf fci + k 2 k\ + k 3 k\ + k 4 

2 6 2fci 1 1 1 

case 2 : — = — — — — , (34) 

/i 2 k k 2 h + k 2 h + k 3 ki + k 4 y ' 

2 _ 4 2fe 4 1 1 1 



10 



A virial theorem proven in [5] asserts that F = 47riV 3 / 2 /i/3, and hence from ^ it 
follows that Vol(y) = 7r 4 /24/^ 2 . We can be more precise about what exactly the seven 
dimensional manifold Y is. It is a tri-Sasaki Einstein space. A tri-Sasaki Einstein space can 
be defined as the level surface of a conical manifold with hyperkahler structure [23]. The 
eight dimensional hyperkahler cone in turn can be thought of as a hyperkahler quotient 
[23] of the the quaternionic manifold H 8 by a subgroup SU{2) x £7(1) 3 C £7(2) x £7(1) 4 . 
We can define the two special £7(1) 's that do not appear in the quotient using the gauge 
connections A a and A of the £7(1) 4 and £7(2) groups respectively. One is the diagonal 
subgroup tr.4+ Yla-A-a- The other is picked out by the CS levels, ktiA + ^2 a k a A a - See 
for example [251 [26] for more details. 

These three cases are related by the Seiberg duality described in the previous section. 
Consider a Seiberg duality on the central node, then on node 1, and then again on the 
central node, SDi SD C SDi. Under such a transformation the CS levels change such that 

(ki,k2,k 3 ,k4) ->■ ^(ki + k2 + k 3 + k 4 ,ki + k2-k 3 -k 4: ,ki-k2 + k 3 -k4 : ,ki-k2-k 3 + k4) . (36) 

Given these relations, it is straightforward to verify that case one is mapped to case two 
and vice versa. The map SD4 SDi SD C SDi takes case three into case one with the rule 

(ki,k 2 ,k 3 ,k4) ->■ ^(ki + k 2 + k 3 + k 4: ,ki + k2-k 3 -k4 i ,ki-k2 + k 3 -k4,ki-k2-k 3 + k l i) . (37) 

To express the action of Seiberg duality in terms of the Weyl group, we need to introduce 
simple roots for D4: a\ = (1,1,0,0), 02 = (—1,1,0,0), a 3 = (0,-1,1,0), and 04 = 
(0,0, —1, 1). Note that a 3 corresponds to the central node. The highest root is then 6 = 
(0, 0, 1, 1) and the CS levels can be expressed in terms of auxiliary variables p a such that 

k\ = pi + P2 , k 2 = -pi +P2 , h = -p 3 + pi , fc 4 = -p 3 - pi . (38) 

The Weyl group acts by permuting the p a and by changing the sign of an even number of 
the p a ^\ We can write a more general expression for fj,, assuming \pi\ > \p2\ > \p 3 \ > \pa\: 

/x 2 2\pi\ + \pi\ + \p 2 \ 2 (\pi\ + \P2\) 2 + ^ 

1 1 



\Pl\ + \P2\ + \P3\ + M M + |P2| + \ps\ - M 
4.1 Comparison with the Numerical Saddle Point 

To check our large N approximation, we can use a computer to calculate the eigenvalue 
distributions numerically. We use the same relaxation algorithm described in [3]. For the 
simple Seiberg dual cases {k\, k2, k 3 , £4) = (2, 0, 0, 0) and (1, —1, —1, —1), we plot the answer 
in figure [2] Note that these cases are somewhat degenerate. Instead of having five distinct 



5 D4 also has an outer automorphism group isomorphic to S3, corresponding to permutations of three 
of the k's. These automorphisms cannot be expressed in terms of Seiberg dualities. Swapping ki and hi 
negates pi; other permutations act less nicely on the pi. Naively, there is an Sa, but the four dimensional 
subgroup generated by the swaps (12)(34), (13)(24), and (14)(23) are Seiberg dualities. 
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Figure 2: The eigenvalue distribution for two Seiberg dual D4 theories with N = 30 
and (fei, ^4) = (2,0,0,0) (left) or (1,-1,-1,-1) (right). The small black points 

correspond to the eigenvalues of the U(2N) node. The large red points correspond to the 
k\ > node. The eigenvalues for the remaining gauge groups - indicated by the blue, green 
and purple points - are coincident. The thin black lines are the large ./V analytic prediction 
with the additional input y « — x/3 (left) and y « x/Q (right). 



piecewise linear regions (for x > 0), we have only two for the (2,0,0,0) theory and one 
for the (1,-1,-1,-1) theory. For (2,0,0,0), the last three regions collapse to the point 
x = /J,/ki, while for (1,-1,-1,-1) the last four regions collapse to the point x = fi/2k\. 

The thin black lines are the analytic, large iV prediction. The agreement is reasonably 
good and should improve as ./V is increased. Note however that we predict only the differ- 
ences between the y a . To make the plot, we first fit y = (yo,i + 2/0,2) /2 to a straight line. 
Given the value of the fit (y « —x/3 for (2, 0, 0, 0) and y x/6 for (1, —1, —1, —1)), the 
transformation (25) is obeyed only up to an overall additive shift of all the y a . 



4.2 Operator counting 

We would like to count the number of gauge invariant chiral operators of the D4 quiver when 
N = 1. To that end, we should first construct words out of the bifundamental fields 
and Bj u that are invariant under the SU{2) C U{2) gauge group. The indices i, j = 1,. . . ,4 
correspond to the four U(l) gauge groups while [A, v = 1, 2 index the t7(2) group. Then we 
quotient by the superpotential (or moment map) relations: 

Y.^B^ + k^ = 0, (40) 

i 

AfB^ + k^i = 0, (41) 
« + = 0. (42) 

We use the convention that we sum over repeated U(2) indices, but we do not sum over re- 
peated U(l) indices. Finally, we restrict to operators with U(l) charges that can be canceled 
by the corresponding diagonal monopole operators which have charge (mki,mk2, rnkj^mk^) 
under the t/(l) 4 , where m is an integer. (There is a sense in which what we do here is a 
generalization of the analysis performed in |27j.) 
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Our counting problem is equivalent to a much simpler counting problem where we ignore 
the (pi and fields. Let R denote the ring of operators. For any .R-module M, let h(M) 
denote the Hilbert series of M. For any ideal /, we have 

h(R) = h(I) + h(R/I). (43) 

If we let / = (</>), then R and / are isomorphic i?-modules and h{((j))) = th{R). So 

h{R) = {l-t)- 1 h{R/{<P)). (44) 

In terms of the function ip(r, m) defined in the introduction, now with a subscript to indicate 
the choice of ring, we have from ^ 

, s fJ- ,(2,1) / rx\ n (ii) / rx\ 

= } { r > 7 J = r ^R/W V 7 "' 7 J • (45) 



When = 0, eqs. (40) and (41) become 

J2 A i B ™ = °> ( 46 ) 

i 

A?B ilt = 0. (47) 

All SU{2) invariant operators may be written as polynomials in AfBj^, e^A^Aj, and 
e^ v Bi^Bj u . These polynomials satisfy various quadratic relations. Some of these relations 
can be obtained by contracting ( |46| ) with bifundamental fields (either Af or e^B^ for upper 
indices, and B^ or c^i/A 1 ? for lower indices). We can obtain other relations by contracting 
bifundamental fields with the identity 

e^£A5 + e^Su + e^e^A = 0. (48) 



Once we apply (47), each of these relations will have zero, two, or three nonzero terms. We 



first mod out by the relations that have two terms. We get two-term relations by contracting 



(46) with two fields with different gauge indices, and by contracting (48) with four fields 
with two or three different gauge indices. 

Modding out by relations with two terms gives us a seven dimensional toric variety. 
The two-term relations are all invariant under a (C*) 7 action. We will let the first four 
charges q±, 92, 93, 94 be the charges under the U(l) gauge groups. We choose 95 so that each 
contraction of an A\ or B\ with an Ai or B2 and each contraction of an ^3 or B3 with an 
A4 or B4 has 95 = 1 and all other contractions have q§ = 0. Likewise, we choose q§ to count 
the number of 1-3 and 2-4 contractions, and 97 to count 1-4 and 2-3 contractions. Note that 
r = 55 + q$ + qj. Also, the following quantitites are even: 

qi + 92 + 93 + 94 , qi + 92 + <?6 + 97 , qi + 93 + 95 + 97 , 91 + 94 + 95 + % ■ (49) 

We can check that all contractions of a pair of ^4's and -B's can be expressed in terms 



of €^ u Bi^Bj v and A^B^v We obtain an expression for A^Bi^ by contracting (48) with two 
A\s, a B2, and a Bi. We can then obtain A^B^ by contracting (46) with Bi and Bk, where 
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l,i,j,k are all different. Finally, contracting (48) with one each of Ai,Bi,Aj,Bj gives us 
e^A^Aj . Therefore, the variety must be seven dimensional. 

For simplicity and to make contact with "case one" discussed above, we will assume 
that 9i > 92 + 93 + 94 , 92> 93 ; 94 > 0. By construction, we have the inequalities 95, 96, 97 > 0. 
We claim the following inequalities also hold, 

95 > 9i + 92 ~ r , 

96 > 9i + 93 - r , (50) 

97 > qi + 94 - r . 

Consider the first of these inequalities. The quantity r is half the total number of A's and 
B's in the operator. The quantity q± + q2 is the number of A\ and A2 fields that we cannot 
pair with with B% and B2 fields respectively. The SU{2) indices of these excess A, fields 
must be contracted in a gauge invariant way. To minimize q§ we can try to contract as many 
of these excess fields as possible with the A3, B3, A4, and -B4 fields. However, if the excess 
number of A\ and A2 fields is larger than r, then we will still have at least 2{q\ + q2 — r) 
A\ and A2 fields left over which we will be forced to contract together, yielding the lower 
bound on q$. 



Now we impose the three-term relations that can be derived from ( 46 ) and ( 48 ) . These 
relations break the symmetries corresponding to 55,(76,97; only the sum q§ + q§ + 97 is 
preserved. We will let 9(x) denote the Heaviside step function, and 9\(x) = x9{x). Any 
monomial that has qj greater than the minimum #1(91 + 94 — r) can be reduced to a sum 
of monomials with q-j = 6\{q\ + 94 — r). 

Therefore, the density of operators with a given q\, 32, 93, 94,?" is 

^J?/(i)( r ' q ^ = \ j dq5 dqe d ( r ~ qi ^ r ~ 95 ~ 96 ~ dl ^ qi + 54 ~ r )) 

• 0(q 5 )6(q 6 )6(q 5 - qi - q 2 + r)9(q 6 - q x - q 3 + r) (51) 

= 2^( r ~ ^)( r ~~ +94~r)- 6»i(gi + q2 - r) - 9 1 (q 1 + 93 - r)). 

The factor of ^ comes from the fact that q\ + q<± + q§ + </6 is always even. If 92 > 93 > 94, 
then in piecewise form we have 



2' 



9i + 92 < r 
^pl, qi + q3<r< qi +q 2 

3r-2q 1 -q 2 -q3 

2 ' 
4r-3gi-g2-93-g4 
2 



9i + 94 < r < qi + 93 (52) 
91 < r < 91 + 94 
0, r < gi 



When we replace qi with ^fcj and multiply by (J,/r, we get agreement with the densities 
computed by the matrix model. 

We may also compare the 5y a j values with numbers of operators in the chiral ring. 
A conjecture in [5] states that for a bifundamental field X ab that transforms in the anti 
fundamental of the a'th U(N) gauge group and the fundamental of the b'th group, 

^2/(]r„ 6 )( r > r!B // i ) = -P( x )[yb( x ) - y a (x) + R(X ab )] (53) 
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where R(X a b) is the R-charge of X a f,. In this case, R(X a b) = 1/2 because of the M = 3 
SUSY. 

First, we count the operators that are not divisible by e^A^A^. Dividing by e^A^A^ 
decreases q$ by one. Since all operators have q§ > 0, any operator with q^ = is not 
divisible by e^A^A^. We can check that all other operators are divisible by e^A^A^. 

Since q$ > q\ + qi — r, operators with q§ = exist only if r > q\ + q2- Therefore, 

^JL^i^) (r ' ?) = 2 9( < r ~ qi ~ q2) ' (54) 

^/(e^MJ) ( r > «) = 2 01 (r " 91 " 92) • (55) 

The factor of 5 comes from the fact that qi + q2 + Q6 + Q7 = Qi + Q2 + r must be even. 
Therefore, operator counting predicts that 

p(x)(l - 2/1 - 2/2 + 2/0,1 + 2/0,2) = (A* - ( k i + fc 2)x). (56) 
We can check that this agrees with our matrix model calculations. Similarly, 

/o(a?)(l -2/1 -2/3 + 2/0,1 +2/0,2) = ^9i(p - (h + k 3 )x) , (57) 
p(x)(l -yi -y 4 + 2/0,1 +2/0,2) = ^i(/w - (fci + A; 4 )x). (58) 

Finally, we count the operators that are not divisible by A^Bip,- Dividing by A^B2^ 
decreases qi, q§ and r by one and increases q2 by one. Since q§ and q§ — q\ — q2 + r must 
both be nonnegative, operators with q$ = or q$ — q\ — q2 + r are not divisible by A^B2^- 

There are two cases to consider. If q§ = 0, then for r > q\ + q2, the value of 
Ygfr^j^Bi ) ( r ' ^) mus ^ be 1/2 as we had before. The second case to consider is q$ = 
qi+Q2-r > 0. This constraint implies that gi+<?2+<?6+<?7 = (<?i+<72-<3 , 5) + (<75+<?6+<77) = 2r 
is even and we do not need an additional factor of 1/2. All operators must satisfy the con- 
straint r > q\. Thus in the range q\ < r < q\ + q2, the value of yr^a^b 2 )( r '^) mus t be 
one. In summary, we find 

^R/(A>tB 2fi )( r > = dl ( r ~ ^l) ~ \ e ^ r ^ ■ ( 59 ) 

Operator counting predicts that 

p(x)(l -yi + y 2 ) =0x(n- hx) - -0 x (n - {h + k 2 )x) , (60) 
and this agrees with our matrix model results. 

5 Discussion 

In this paper, we saw that there is an ADE classification of M = 3 SUSY CS theories 
in 2+1 dimensions with bifundamental matter fields whose partition function Zgz has a 
nice large N limit. By nice limit, we mean that when we use the matrix model of f2] 
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Figure 3: We conjecture that these are the six families of quivers involving both unitary 
and orthogonal/symplectic gauge groups that have M = 3 supersymmetry. The circular 
nodes correspond to unitary groups, and the square nodes correspond to orthogonal or 
symplectic groups. These quivers can be interpreted as non-simply-laced extended Dynkin 
diagrams. An edge connecting a unitary gauge group to an orthogonal/symplectic gauge 
group corresponds to a double bond in the Dynkin diagram, with the arrow pointing toward 
the unitary gauge group. We are not sure of the precise rules for deciding which gauge groups 
should be orthogonal and which should be symplectic. There does not appear to be a way 
to realize the extended Dynkin diagrams G2, Ii, ^2^-1 > or ^4 as <l u i vers > since we do not 
have a way of interpreting triple or 00 bonds, or a pair of double bonds pointed in the same 
direction. 

to calculate Zg3, there is cancellation in the long range forces between the eigenvalues. 
The classification raises a number of questions which we were not able to answer in this 
paper. Leading among them is the significance of the Nj = 2N C condition Q. For a 3+1 
dimensional gauge theory, this condition implies the vanishing of the NSVZ beta function 
and the existence of a conformal fixed point, but in the present context the meaning is 
obscure. It is tempting to speculate that Nf = 2N C is necessary in 2+1 dimensions for 
the existence of the fixed point. However, it could be that the condition is just a technical 
limitation of the method introduced in [3] for calculating the large N limit of the matrix 
model. 

Another question has to do with the restriction to U(N) gauge groups. It would be 
interesting to study quivers of SO(N) and Sp(N) groups in the large N limit. We conjecture 
that "good" theories can all be described by extended Dynkin diagrams. We believe that the 
"good" quivers containing just SO(N) and Sp(N) gauge groups are simply laced extended 
Dynkin diagrams, while the "good" quivers that also contain U(N) are the ones shown in 
Figure [3j It would be interesting to study this correspondence further. We note that an 
SO-Sp version of the simply laced A\ quiver was studied in [28J . 
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We would like to end this paper by describing some progress we made in writing F in 
a way that is manifestly invariant under Seiberg duality. In the introduction, we reviewed 
the fact that Vol(Y) ~ 1/F 2 where Y is the tri-Sasaki Einstein space in the dual eleven 
dimensional supergravity description of the CS theory. In the A n and D4 cases, we have 
found that Vol(Y) can be written as a rational function of the CS levels k a . The numerator 
of this function takes the form 

n 

^2 det(«i, . . . ,a n ) 2 Y[ \a a ■ p\ (61) 

Ql,...,Q n G-R a=l 

where R is the set of roots of the A n , D n , or E n Lie algebra. Recall that the CS levels can 
be reconstructed from the simple roots via k a — OL a • p. 
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A A Simple Class of Solutions 



In this section, we find a simple class of solutions for any ADE quiver. Let us start by 
looking at the Eq quiver. 

At the risk of unacceptably proliferating the number of indices, we label the CS levels of 
the external nodes k\ a , of the internal nodes k 2a , and of the central node k^ with the con- 
straint J2 a (kia + 2k 2a ) = — 3&3. The imaginary parts of the eigenvalues are correspondingly 
labeled yi a , y 2a ,i and y^,j- For small enough x, the equations of motion yield 



p(x) 



ft 
6 



Via ~ V2a,I 



3k\ a x 



U2a,I - U3,J 



(k la + 2k 2a )x 



(62) 



By small enough, we mean that no difference between the y's has saturated to ±1/2. More 
precisely, we require \x\ < n/6\ki a \ and \x\ < fi/2\k\ a + 2/c2 a | for all a. 

As in the D4 case, for a certain class of solutions the eigenvalue distribution ends inside 
this central region. Let us try to extremize the free energy as a function of the endpoint 
. The on shell free energy as a function of x* is 



F(x*) = N 3 / 2 
which leads to the condition that 



(63) 



J2(Zk 2 la + (k la + 2k 2a ) 2 ) 



-1/4 



(64) 
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In order for to lie in the central region, we must take k\ a — k^b — ±& or e 2k\ a — — kib — 
±2k, leading to the final result 

F = TrN 3/2 A^/2\k\ , (65) 

and the result for the volume 

Vol < F > = 4^ ' (66 > 
The tri-Sasaki Einstein manifold in question is a level surface of the hyperkahler quotient 
M M ///SU(3) x SU(2) 3 x U(lf. 

The solution above has some characteristic features that we can try to generalize for 
arbitrary ADE quivers. The y\ a — y2a,l and y2 a ,l — U3, J ah saturate to ±1/2 at x = x*. 
Also, the absolute values of the 5y are equal, \y\ a — V2a,i\ = \V2b,i — U3,j\- For the arbitrary 
ADE case, let the CS level of the U(n a N) gauge group be k a = sgn(a)kn a , where we choose 
sgn(a) = ±1 so that neighboring gauge groups have opposite sign. This choice of CS levels 
is consistent with the condition ]P a n a k a = 0. The equations of motion following from 



(21) are satisfied if we set all the \Sy\ equal to each other and the eigenvalue density to a 



constant: 

Sy a b,u = sgn(a)— , p = = — 2 • ( 67 ) 

With a little more work, we can find the chemical potential and the volume of the associated 
tri-Sasaki Einstein manifold 
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